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Abstract. In the context of nanowire heterostructures we perform a discrete 
to continuum limit of the corresponding free energy by means of P-convergence 
techniques. Nearest neighbours are identified by employing the notions of 
Voronoi diagrams and Delaunay triangulations. The scaling of the nanowire is 
done in such a way that we perform not only a continuum limit but a dimension 
reduction simultaneously. The main part of the proof is a discrete geometric 
rigidity result that we announced in an earlier work and show here in detail for 
a variety of three-dimensional lattices. We perform the passage from discrete 
to continuum twice: once for a system that compensates a lattice mismatch 
between two parts of the heterogeneous nanowire without defects and once 
for a system that creates dislocations. It turns out that we can verify the 
experimentally observed fact that the nanowires show dislocations when the 
radius of the specimen is large. 


Introduction 

Rigidity results of elastic materials have been of great interest in mathemati¬ 
cal continuum mechanics in recent years, in particular since the seminal work by 
Friesecke, James and Muller [6]. Such results yield a deeper insight into the prop¬ 
erties of materials through an estimate of the distance of the deformation gradient 
from the set of rotations; this distance is in turn estimated from above by the free 
energy of the system. 

The rigidity estimates turn out to be crucial steps in various proofs as for instance 
of P-convergence results in the context of dimension reduction. This was also the 
case in our earlier paper [11], in which we derived a discrete to continuum limit 
and a dimension reduction of an energy of a heterogeneous nanowire (see [12] for 
an abridged version). There we presented a detailed analysis of the passage from 
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the two-dimensional setting to the one-dimensional limit, and we gave a summary 
of the corresponding dimension reduction from three dimensions to one dimension. 
The purpose of this article is to show the rigidity estimates (Section 3) and the 
main features of the latter case in detail. 

Further, we elaborate on various three-dimensional lattices that are of impor¬ 
tance in applications: the face-centred cubic lattice, the hexagonal close-packed, 
the body-centred cubic lattice and the diamond cubic lattice, see Section 1. These 
lattices occur for instance in aluminum and gold, magnesium and zinc, iron and 
tungsten, and germanium and silicon, respectively. Note that Si/Ge nanowires 
have applications in the semiconductor optoelectronics [10, 15]. In Section 3 we 
show that our discrete rigidity result applies to all these lattice structures. The 
main property of such lattices is their geometric rigidity: they define a tessellation 
of the space into rigid polyhedra whose edges correspond to bonds in the lattice. 
Our approach does not work in non-rigid lattices, like a simple cubic crystal with 
nearest-neighbour interactions only. 

We are interested in the mathematical modeling of dislocations in heterogeneous 
nanowires. We assume that the material consists of two parts with the same lattice 
structure but different lattice constants. The interface between the two parts is 
assumed to be flat. The material overcomes the lattice mismatch either defect-free 
or by creating dislocations. As was pointed out by Ertekin et al. [3], it is the radius 
of the nanowire which determines whether the material creates dislocations or is 
defect free. In our model the radius roughly corresponds to the number of layers 
of atoms parallel to the direction of the wire, see Section 2. We prove that it is 
energetically more favourite to create dislocations than to relieve the mismatch in a 
defect-free way if the thickness of the nanowire is sufficiently large (see Remark 7). 

The underlying idea of our mathematical model, which we introduce in Section 2 
in detail, goes back to the variational model proposed in [13] in the context of non¬ 
linear elasticity and which was later generalized to a discrete to continuum setting 
in [11]. As before we assume that the total energy only consists of nearest-neighbour 
interactions which are harmonic, though it is possible to generalize this as discussed 
in [11, Section 4]. In order to be able to apply a rigidity estimate, we always im¬ 
pose a non-interpenetration condition, which ensures that the deformations of the 
discrete setting preserve the orientation of each cell; similar assumptions were made 
e.g. in [2, 7]. The non-interpenetration assumption can be dropped if one takes into 
account interactions beyond nearest neighbours, see the recent work [1]. It is worth 
mentioning that a related variational model for misfit dislocations has been recently 
proposed in [4]. 

As in [11] we distinguish the systems with and without defects already in the 
given reference configuration. For both such systems we study the corresponding 
free energy of nearest- neighbour interactions in a discrete to continuum limit with 
dimension reduction. For the definition of the nearest neighbours in the discrete 
settings close to the interface it is useful to work with the notion of Delaunay 
triangulations and Voronoi cells, see [11], where this was introduced for the first time 
to describe configurations with dislocations, see also Section 1 for an introduction. 

In Section 4 we compare the minimizers of the limiting functionals, which char¬ 
acterize the minimum cost needed to compensate the lattice mismatch with and 
without defects, respectively. It turns out that this cost depends on the thickness of 
the wire described by a mesoscale parameter k. More precisely, it depends quadrat- 
ically on k if there are dislocations, and scales like k^ if there are no defects. Hence 
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for sufficiently large k, i.e. large radius of the wire, dislocations are energetically 
preferred. The result is based on a scaling argument. In particular for applications 
in semiconductor optoelectronics it would be interesting to know the threshold kc 
below which the nanowire deforms defect-free. This is however out of reach with 
our current methods so that we leave this as an open problem for future research. 

1. Three-dimensional lattices 

We consider various three dimensional lattices whose unit cells are rigid convex 
polyhedra. In this context, rigidity is understood in the following sense: once the 
lengths of the edges of a polyhedron are given, then the polyhedron is determined 
up to rotations and translations, under the assumption that the polyhedron itself 
is convex. We recall that a convex polyhedron is rigid if and only if its facets are 
triangles, according to the classical Cauchy Rigidity Theorem (see, e.g., [9]). We 
consider four types of discrete lattices in dimension three: the face-centred cubic, 
the hexagonal close-packed, the body-centred cubic, and the diamond cubic. They 
should be interpreted as prototypes to which our approach can be applied, under 
slight modifications in each case. For a general overview on lattice structures see, 
e.g., [8]. 

All the lattices we will introduce, fulfil a property of rigidity. Indeed, the cor¬ 
responding nearest-neighbour bonds provide a tessellation of the space into rigid 
convex polyhedra, as we will make precise case by case. (In the diamond cubic, also 
next-to-nearest neighbours will be used.) We always assume a non-interpenetration 
condition, see (7) below. 

A major role in modeling is then played by the choice of the nearest neighbours of 
each lattice. Here they are defined according to the notion of Delaunay pretriangu¬ 
lation, as given in the following definitions. Such a general definition can be applied 
also when the lattice is irregular, so in particular across the interface between the 
phases, see Section 2.1. 

For later convenience we give the definition for all dimensions N > 2. Let C C 
be a countable set of points such that there exist R,r > 0 with inf^-ggAr ^(£ fl 
B{x,R)) > 1 and \x — y\ > r for every x,y G £, x ^ y, where B{x,R) := {y G 
|a;- 2 /| < R}. 

Definition 1.1 (Voronoi cells). The Voronoi cell of a point x G C is the set 

C{x) := {z G : \z — x\ < \z — y\ Wy G £} . 

The Voronoi diagram associated with £ is the partition {C{x)}x^c- 

Definition 1.2 (Delaunay pretriangulation). The Delaunay pretriangulation asso¬ 
ciated with £ is a partition of in open nonempty hyperpolyhedra with vertices 
in £, such that two points x,y G £ are vertices of the same hyperpolyhedra if and 
only if C{x) fl C{y) ^ 0. 

Definition 1.3 (Nearest neighbours). Two points x,y G £, x ^ y, are said to be 
nearest neighbours (and we write: x, y NN) if they are vertices of an edge of one of 
the hyperpolyhedra of the Delaunay pretriangulation. 

Definition 1.4 (Next-to-nearest neighbours). Two points x,y G £, x ^ y, are said 
to be next-to-nearest neighbours (and we write: x, y NNN) if, setting 

£^{x) := £\{y: X, y NN} , 
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we find that fl C^{y)) > 0, where {C^{y)}y£C,(x) is the Voronoi dia¬ 

gram associated with £* (x). 

The Voronoi diagram and the Delaunay pretriangulation associated with a lattice 
are unique. For these and other properties we refer to [11, Section 1] and references 
therein. 


1.1. FCC lattice. The face-centred cubic lattice is the typical structure of metals 
such as aluminium, gold, nickel, and platinum. It is the Bravais lattice generated 
by the vectors 

vf :=v^(l,0,0), vf :=y2(i,i,0), vf := ^2(0, i, i) , 

namely 

■= {fivf +6^2 : ^ 1 ,^ 2 , ^3 S Z}. 

The resulting lattice is obtained by repeating periodically in the space a cubic cell 
of side where the atoms lie at the vertices and at the centre of each facet. It 
is readily seen that two points x,y G are nearest neighbours in the sense of 
Definition 1.3 if and only if \x — y\ = 1, i.e., they are joined by half a diagonal of 
a facet of the cubic cell. Each atom has twelve nearest neighbours. The Delaunay 
pretriangulation provides a subdivision of the space into regular tetrahedra and 
octahedra of side one, thus in rigid convex polyhedra, see Figure 1. Remark that 
the diagonals of the octahedra, whose length is v^, correspond to next-to-nearest 
neighbours. The latter will not enter the definition of the energy (3). 


1.2. HCP lattice. Our approach works also for non-Bravais lattices such as the 
hexagonal close-packed structure found in some metals as, e.g., magnesium and 
zinc. It is defined by 

^ '■= {uf + Civf -I- + ^3^3 : ^1, ^2: ^3 € Z, i = 1, 2} , 


where 



are generators of two sublattices and 


1 VI 

2 ’ 2 


0), 



1 

2 ’ 2 


0 ) 


uH:= (0,0,0), uH:=(0,^,^) 

are called vectors of the basis. The lattice is thus obtained by merging two Bravais 
sublattices (defined for i = 1 and i = 2, respectively). As in the previous case, the 
nearest neighbours are those couples with distance one, each atom has twelve nearest 
neighbours, and the Delaunay pretriangulation consists of regular tetrahedra and 
octahedra of side one, see Figure 2. As before, the diagonals of the octahedra, which 
correspond to next-to-nearest neighbour interactions, will not enter the definition 
of the energy (4). 


1.3. BCC lattice. The body-centred cubic lattice is typical of some metals as, 
e.g., iron and tungsten. It is the Bravais lattice generated by the vectors 


/f := ^(-1,1,1) 


._ VI 


:=^(1,-1,1), vf :=^(1,1,-1) 


namely 

■= {^ivf -I- ^ 2^2 + Csvf : ^ 1 ,^ 2 , ^3 S Z} . 

The resulting lattice can be viewed by repeating periodically in the space a cubic 
cell of side where the atoms he at the vertices and at the centre of the cube. 
According to Definition 1.3, the nearest neighbours are those couples with distance 
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Figure 1. In the face-centred cubic lattice the nearest-neighbour 
structure of the atoms provides a subdivision of the space into 
tetrahedra (a) and octahedra (b). Figure (c) shows a quarter of an 
octahedron in the same unit cell. Grey dots denote points lying on 
the hidden facets. 




Figure 2. The hexagonal close-packed lattice is associated with 
a tessellation of tetrahedra and octahedra as the ones in the fig¬ 
ure. Only some of the bonds and some of the polyhedra of the 
pretriangulation are displayed. 



Figure 3. The body-centred cubic lattice is associated with a tes¬ 
sellation of irregular tetrahedra as the one in the figure. 
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(i.e., those joined by half a diagonal of the cubic cell,) as well as those couples 
with distance (i.e., those joined by an edge of the cubic cell). Thus, in contrast 
with the face-centred cubic, in this case the notion of nearest neighbours differs from 
other notions based on the Euclidean distance. According to this definition each 
atom has 14 nearest neighbours. Correspondingly, the Delaunay pretriangulation 
consists in a subdivision of the space into irregular tetrahedra, with four edges of 
length and two of length see Figure 3. Such an asymmetry in the definition 
of nearest neighbours leads to consider an anisotropic energy, see (5). 

1.4. DC lattice. Finally, we present the diamond cubic lattice, which is composed 
of two interpenetrating face-centred cubic lattices (thus, it is non-Bravais). It is 
relevant in applications to nanowires, since it is the structure of materials of use, 
such as silicon and germanium [10]. When the sites of the two interpenetrating lat¬ 
tices are filled with two different species of atoms, the structure is called zincblende 
and is typical of Gallium arsenide (GaAs) and Indium arsenide (InAs), also used in 
technical applications to semiconductor optoelectronics [10]. 

The diamond cubic structure is defined by 

CP := {uf -f Civ° -f C 2 V° + 6,6,6 S Z, i = 1,2}, 

where v^ := vj, j = 1, 2,3, are as in the face-centred cubic and 
u? := (0,0,0), u? := 72(1,1,1) 
compose the basis. It is convenient to split the lattice as follows, 

( 1 ) 

:= {uf -h 6v° -h 6v° -H 6v° : 6,6,6 S Zj , i = 1, 2 , 
where the sublattices i = 1,2, are face-centred cubic, see Figure 4. 

Each atom of the sublattice x G £^',i = 1,2, has four nearest neighbours at 
distance all belonging to the sublattice £^^ ,j ^i. Such bonds are not enough to 
provide a rigid tessellation of the space. Therefore we need to take into account also 
the next-to-nearest neighbours. By Definition 1.4, the next-to-nearest neighbours 
of X in £P turn out to be its nearest neighbours as an element of More 

precisely, each atom x lies at the barycentre of a tetrahedron whose vertices are 
the nearest neighbours of x; the edges of such a tetrahedron correspond to next-to- 
nearest bonds. Thus, when next-to-nearest neighbours are considered, £^ inherits 
some rigid structure from the (face-centred cubic) sublattices i = 1, 2. 

For a better understanding of the diamond cubic lattice, we also refer to the 
simpler example of the planar honeycomb lattice, which can be treated by the 
same methods as presented here. This two-dimensional example contains the main 
ideas for treating non-Bravais lattices with next and next-to-nearest neighbours, see 
Figure 5. 


2. Setting of the model 

In order to mathematically describe the three-dimensional heterostructured nanowires 
we introduce four parameters e, k, X and p, next to the lattice structures discussed 
above. 

The parameter e > 0 scales the equilibrium lattice distances and allows con¬ 
sidering a passage from the discrete to the continuous setting by letting £ —>• 0+. 

The parameter k € N, k > 1, mimics the thickness of the nanowire. The shape 


RIGIDITY OF THREE-DIMENSIONAL LATTICES 


7 



Figure 4. Cubic cell in the diamond lattice Atoms from the 
sublattice are represented in black/grey, while white atoms are 
from the sublattice Nearest-neighbour bonds are displayed 

by solid thick lines. Moreover, the picture shows a tetrahedron 
from the Delaunay pretriangulation of its edges (solid and 
dashed thin lines) correspond to next-to-nearest neighbours in . 
A white atom lies at the barycentre of the tetrahedron, which is 
further divided into four irregular tetrahedra by the bonds between 
the barycentre and each vertex. 



Figure 5. Bonds and triangulation in a honeycomb lattice. The 
lattice is given by := U where := {uf -|- ^iv* -|- 
6,6 e Z}, v? := (1,0), v| := (i,^), u? := (0,0), 
u* := (0, ■^). This results into two interpenetrating sublattices 
and both being hexagonal (i.e., equilateral triangular). Atoms 
from and are displayed in different colors in the picture, 
respectively in black and in white. In the left part of the figure we 
indicate nearest neighbour (solid) and next-to-nearest neighbour 
bonds (dashed lines). The right part of the figure shows a possi¬ 
ble triangulation, that is the natural triangulation of enriched 
by considering the nearest-neighbour bonds between atoms x G Cf 
and y G . This corresponds to ignoring the bonds between atoms 
of £2 , cf. Section 2.3. 
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of the nanowire in the discrete setting is a parallelepiped of length 2 L, L > 0 , 
and width and the height fee, see Section 2.2 for details. In the continuum limit 
e —O'*', the length is conserved whereas the width and height tend to zero thus 
giving a dimension reduction of the system from three to one dimension. Still, the 
microscopic parameter k has an impact on the continuum energy, which then allows 
investigating the limiting behaviour in dependence of the microscopic thickness k 
of the wire. 

The parameters A and p allow modeling the microscopic biphase structure of the 
nanowire. Here, A € (0,1) denotes the ratio of the equilibrium distances in the 
deformed configuration of the material on the right hand side of the interface and 
of the material on the left hand side of the interface, see Section 2.1 for details. 

The parameter p G (0,1] gives the ratio of the lattice distances of the two parts of 
the material in the reference configuration, where p G [A, 1] is the most interesting 
case. This allows treating different geometries of the nearest neighbours and in 
particular for dislocations. The case of a defect-free body is modeled by p = 1; 
the coordination number, i.e., the number of nearest neighbours of any internal 
atom, is constant in the lattice. If the crystal contains dislocations in the reference 
configuration, the coordination number is not constant. As we will show, this is the 
case for p 1 and k sufficiently large. 

2.1. Biphase lattices and rigid tessellations. Given p G (0,1] and vectors 
vi,V 2 ,V 3 , ui,U 2 G we define the biphase atomistic lattice 

( 2 ) £,:=£-U£+ 

by juxtaposing the two lattices and £+ given by 

£i := {ui -I- ^ivi -I- ^2V2 + C 3 V 3 : ^ 1 , ^ 2 , ^3 € Z , i = 1, 2 , < 0} , 

•= + ^I’^l + ^2^2 + ^3^3 : ^1, ^3 € pZ , f = 1, 2 , ^1 > 0} . 

We will apply the above definitions to the crystals introduced in Section 1 and denote 
by £p, £p , £p, and £^ the lattices obtained by taking the vectors {vj, uf }j=i, 2 , 3 , 

{v? uf}j=i, 2 , 3 , {v? uf }j=i, 2 , 3 , and {v°, uP}j=i, 2 , 3 , respectively, where uf = uf := 

i=l,2 i=l,2 i=l,2 

0. 

In each of the four cases we find similar structures for the planes at the interface 
between the lattices and £+. More precisely, for the face-centred cubic and 
the hexagonal close-packed, the interfacial planes are two-dimensional equilateral 
triangular Bravais lattices, see Figure 6. In the body-centred cubic, the interfa¬ 
cial planes are triangular Bravais lattices, but not equilateral, since the distance 
between nearest neighbours is not constant. Finally, in the diamond cubic, whose 
properties are similar to the face-centred cubic, we also find equilateral triangular 
planes composed by atoms of one of the sublattices. (For a lower dimensional idea, 
see Figure 7.) 
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Figure 6. By cutting a cubic lattice along certain transverse 
planes, one finds two-dimensional hexagonal Bravais lattices, (a) 
face-centred; (b) body-centred. 



Figure 7. Dislocations in a honeycomb-type lattice. The bonds 
at the interface are chosen in the following way: First one consid¬ 
ers only black atoms and finds a Delaunay pretriangulation, which 
is then refined to a triangulation (dashed lines); the same is done 
for white atoms (dotted lines). The dashed and dotted lines thus 
obtained give the bonds between next-to-nearest neighbours. Fi¬ 
nally, each white (resp. black) atom lying inside a triangle formed 
by three black (resp. white) atoms is connected to the vertices of 
that triangle by nearest-neighbour bonds (solid lines). 
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Next we define the interfacial bonds in the case when p ^ 1- Following the idea 
already used in the regular parts of the lattices, we consider the (unique) Delaunay 
pretriangulation Tp of Lp (Definition 1.2). This defines, in the case of and 

, a tessellation of the space into rigid polyhedra away from the interface. At the 
interface, the partition Tp may contain polyhedra with quadrilateral facets (to see 
this, one should recall that the interfacial atoms lie on two parallel planes consisting 
of two-dimensional triangular Bravais lattices, with parallel primitive vectors): in 
such a case we refine Tp further, in order to obtain rigid polyhedra. More precisely, 
given a quadrilateral facet we introduce a further bond along a diagonal of the 
facet; correspondingly, the region around the interface is subdivided into (irregular) 
tetrahedra and octahedra. 

Following this construction we define a partition of the space into rigid polyhedra 
and call it the rigid Delaunay tessellation associated to Cp, denoted by Tp. The 
nearest neighbours are the extrema of the edges of the polyhedra of the subdivision. 
Such procedure can be followed for and . Instead, in the diamond 

cubic lattice, applying Definition 1.3 may result in nearest-neighbour bonds between 
interfacial atoms of the same sublattice (which should instead be next-to-nearest 
neighbours). This would not be consistent with the structure defined away from 
the interface, therefore we follow a different construction. 

Recall that consists of two interpenetrating face-centred cubic lattices 
and see (1) and (2). We introduce Delaunay pretriangulations for the sublat¬ 
tices ^ and individually, which we further refine in order to obtain triangular 
facets as before; we say that the vertices of the resulting edges are next-to-nearest 
neighbours in . The tessellation of consists of (possibly irregular) tetrahedra 
and octahedra; some of them may contain one atom x of the other sublattice £p ^, 
j ^ i. In this case, we connect x to the vertices of the surrounding polyhedron 
and say that each of those vertices is a nearest neighbour for x. When applied to 
the regular parts of the lattice, this construction is consistent with the notion of 
nearest and next-to-nearest neighbours presented in the previous section. (For a 
simpler idea about the resulting structure, we refer to Figure 7 in the case of a 
honeycomb-type lattice.) 

2.2. Reference configurations and interaction energies. We now pass to 
rescaled, bounded lattices. Given L > 0, e G (0, 1], and fc G N, we define 

£p,e(A:) := (e£p) n , 

where 

:= {^ivi + ^2V2 -I- ^3^3: Cl S (— L, L) , C 2 , Ca S (0, ke)} 

and DfeE is the union of all (closed) polyhedra of Tp^e that intersect D.ke on a 
nonempty set (see (6) for the definition of Tp^e)- Notice that the lattices £p£, 
£p_£, and denote the corresponding rescaled, bounded lattices for the 

other crystal structures, with the lattice vectors chosen as above. We set 

•= {Civi + C 2 V 2 + C 3 V 3 : Cl S (-£,0), C2,C3 G (0, ke)} , 

^ke ■= {Civi + C 2 V 2 + C 3 V 3 : Cl S (0, L) , C 2 , Ca S (0, ke)} , 

■= {eu* -b CiVi + C 2 V 2 + C 3 V 3 G £p,e{k ): Cl < 0} , 

:= -b Civi -b C 2 V 2 -b C 3 V 3 G Cp,e{k ): Cl > 0} . 
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Two points x,y G Cp^e are said to be nearest (resp., next-to-nearest) neighbours if 
a;/e, y/e fulfil the corresponding property in the lattice Lp. This definition applies 
to each of the four cases presented above. Notice that we denote these lattices by 

and as well as by 




and 


Cp^^ (fc), respectively. 

We have introduced so far the bonds that enter the definition of the energy, 
which we generally denote by . Next we specialise £^ for each of the four lattices 
introduced above. In the cases of the face-centred cubic and of the hexagonal close- 
packed, the total interaction energy is defined respectively by 


£j'^{ue,p,k) := 


(3) 


^ E 

X, y NN 

yecl^ik) 


Ue{x) - Ue{y) 


- 1 


+ h 


E 

X, y NN 


Ue{x) - Ue{y) 


-A 


ye£ 


for every deformation u,,: 


and by 


7H, A 


ius,p,k) := 


(4) 


^ E 

X, y NN 
(fe) 

+ ^ E 

X, y NN 

(k) 

vecl,{k) 


Ue{x) - Ue{y) 


- 1 


Ueix) - Ue{y) 


-A 


for every deformation u^: C^^g{k) —>• R^. 

For the body-centred cubic, we need to use an anisotropic energy, because of the 
different length of the bonds between nearest neighbours in the reference configu¬ 
ration. For every deformation —>■ R^ we define 


(5) 


^B,A 


{Ue,P,k) 


\ E 

X, y NN 

xG£?,7(fc) 


Ue{x) - Ue{y) 

e 



2 





where ipy := ip{ |^_^| ) and c/?: —>■ (0, +cx)) is a smooth function snch that 

V5(v) = V2 if V e {(±1,0, 0), (0, ±1, 0), (0, 0, ±1)} , 

(^(v) = ^ if Vsve {(±1,1,1),(±1,-1,1),(±1,1,-1),(±1,-1,-!)}. 
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Finally, recall that the diamond cubic lattice consists of two interpenetrating 
face-centred cubic lattices. Therefore we set for a deformation Ue' C^g{k) — 


7D. A 


{ug, p, k) :— -f C2 £'p^^j,j + £] 


"NN J 


where the first two summands account for next-to-nearest neighbour interactions 
and are defined as in (3), namely 


^NNN 


i E 

X, y NNN 
(fc) 



D . 

x^Cp'e (fc) 


while the last term is 

Ue{x) - Ue{y) 


cD 

'NN 


E 

X, y NN 
(^) 

yec^Ak) 


4 


i E 

X, y NN 

(k) 


Ue{x) - Ueiy) 


\ V6 
- ^ — 


The choice of the constants ci, C 2 > 0 determines how strong the interactions be¬ 
tween atoms of the same sublattice are. 


2.3. Admissible configurations. In order to define the admissible deformations, 
we introduce piecewise affine functions. To this end, we need to refine 7^ to a proper 
triangulation. However, we do not change the definition of the nearest neighbours, 
i.e., we do not introduce new interactions in the energy. 

Remark 1. For the reader’s convenience, we summarise here the different tes¬ 
sellations of the space associated to a biphase discrete lattice £p, adopted in our 
setting. 

• We have started from the (unique) Delaunay pretriangulation Dp (Definition 
1.2), which may contain non-rigid polyhedra at the interface. 

• We have refined Dp, obtaining a rigid Delaunay tessellation Dp, a partition of 
the space into (possibly irregular) tetrahedra and octahedra. Such a tessel¬ 
lation is not unique, indeed we have chosen a diagonal for each quadrilateral 
facet of polyhedra of Dp- The corresponding bonds enter the definition of 
the interaction energy. 

• In order to work with piecewise affine functions, in this section we further 
refine Dp to get three possible triangulations (i.e., subdivisions of the space 
into tetrahedra only), denoted by Dp'^\ and D}'^\ respectively. 

The above construction is used to work in the case of C^, and For £^, the 
definition of Dp is different, as made precise in Sections 2.1 and 2.3. 

In the case of £p, £^, and £^, given a (possibly irregular) octahedron of 7^, we 
divide it into four irregular tetrahedra by cutting it along one of the three diagonals. 
We choose the diagonal starting from the vertex with the largest xi-coordinate; if 
two or three vertices have the same largest xi-coordinate, we take among them the 
point with largest X 2 -coordinate; if two of such vertices have also the same largest 
a; 2 -coordinate, we take the one with the largest Xa-coordinate. By repeating the 
process on every octahedron of Dp, we obtain a triangulation that we denote by 
Dy^. Other two triangulations 7)^^^ and D}^^ are obtained by repeating the same 
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procedure, but with different ordering of the indices, namely X 2 ,X 3 ,xi and X 3 ,xi,X 2 
respectively. 

In the case of the diamond-cubic lattice, we define a triangulation as follows: we 
consider the Delaunay pretriangulation of , which is rigid. As already observed, 
some of the tetrahedra of the latter pretriangulation contain an atom of at the 
barycentre (more precisely, every other tetrahedron has this property, see Figure 
4). Such tetrahedra are further subdivided by connecting the barycentre to the 
vertices. In other words, we define a tessellation into tetrahedra and octahedra 
by considering the (nearest neighbour) interactions between atoms x G and 
y € as well as the interactions between atoms of (nearest neighbour 

if restricted to next-to-nearest neighbour if viewed in the whole £°), and 

ignoring the interactions between atoms of We apply the same rule to the 

biphase lattice and further subdivide the resulting octahedra as done for £p, 

, and £^, obtaining three possible triangulations. For a better understanding we 
illustrate the tessellation thus defined in the simpler case of the honeycomb lattice 
in Figure 5. 

Given a function u: £p —>■ , we denote by and its piecewise affine 

interpolations with respect to the triangulations 'Tp'^\ find 7'p'^\ respectively. 

Analogous definitions and notations hold for the rescaled bounded lattices of the 
type £p^e{k). More precisely, define 

( 6 ) Tp,e--={eT-.TGrp} and TI;$ := {eT: T 

for i = 1, 2, 3. The set of admissible deformations is 

itg piecewise affine, 

Vue constant on Vtkg fl T V T G Tp^J , 
det Vug > 0 a.e. in , 

Ue(V) is convex WV G 7^,e} • 

The restriction of Ue G Ap^e(^ke) to £p^e{k) is still denoted by u^- We will see that 
the limiting functional is independent of the choice of the triangulation I'p^J in (7), 
cf. Remark 5. 

Remark 2. The assumption of convexity on the images of the octahedra of Tp^e 
is needed to enforce rigidity: without such an assumption an octahedron could be 
compressed without paying any energy. On the other hand, the notion of non¬ 
interpenetration used in (7) is independent of the choice of the triangulation 
provided the image of each octahedron is assumed to be convex, as clarified by 
Lemma 3.3. 

It will be convenient to introduce 

:= { 6 vi + 6 v 2 +C 3 V 3 : Cl S (-oo,-hoo), ^ 2,6 e ( 0 , k)} 

and to denote by flk,oo the union of all (closed) polyhedra of Tp that have a nonempty 
intersection with £tk,oo- We define the set of admissible deformations on the rescaled 


Ap,e{^ke) :={u. GG°(Ofc,;K3): 

(7) 
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infinite domain as 

-^p.oo(f^fc.oo) := {m G C'°(r2fc^oo; : u piecewise affine, 

Vu constant on nk,oo fi T VT G 7^ , 
det Vm > 0 a.e. in f2fc_oo , 
u{V) is convex VP G 7^} . 

All definitions apply to each of the four cases presented above. Correspondingly, we 
define the energy on the rescaled infinite domain and denote it by . Specifically, 
given a discrete deformation v of the face-centred cubic lattice, is defined by 

(8) 

:=\ ^ (|w(x)-i;(2/)| - 1)^-h i ^ {\v{x) - viy)\ - X)"^ . 

X, y NN X, y NN 

x€Cf {k) xGC^'^ (k) 

yeC^ik) 

x,y&nk,^ x.yGfifc.oo 

Analogous definitions hold for C^, and 

3. Discrete rigidity in dimension three 

A key tool in the analysis developed in [11] for two-dimensional heterogeneous 
nanowires as well as in the analysis of the three-dimensional setting is the following 
rigidity estimate. 

Theorem 3.1. [6, Theorem 3.1] Let N > 2, and let 1 < p < -|-oo. Suppose that 
U C is a hounded Lipschitz domain. Then there exists a constant C = C{U) 
such that for each u G there exists a constant matrix R G SO{N) such 

that 

(9) 11 Vm - Rh^UM^x^ < C(t/)lldist(Vw, 50(iV))II. 

The constant C{U) is invariant under dilation and translation of the domain. 

In order to employ the above result, we need the discrete rigidity estimates of 
Lemmas 3.2 and 3.4, which state that the energy of a lattice cell is bounded from 
below by the distance of the deformation gradient from the set of rotations. Similar 
rigidity estimates are used in [2, 5, 14, 16]. 

We use the following notation for the vectors determined by the edges of the 
regular tetrahedron S of edge length one: wi := (1,0,0), W2 := (|,-^,0), W3 := 
W2 — wi, W4 := (i, ■^, ■^), W5 := W4 — W2, and we := W4 — wi, cf. Figure 8. 

Lemma 3.2. There exists C > 0 such that 

6 

(10) dist2(F,5'0(3)) < C^(|Fw,| - 1)2 

i=l 

for every F G GL'^{3). 

Proof Set := |Pwi| -1 and S := (Si,..., Se), then “ 1)^ = SLi = 

|(5|^. Without loss of generality we may assume that 

Fwi = (1 -I- i5i)wi, Fw2 G span{wi,W2} , Fw 2 ■ 62 > 0 , 
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Figure 8. The tetrahedron S and its image F{S). 


as in Figure 8 . Notice that the above assumptions imply Fw 4 -63 >0. We have 
dist2(F, 50(3)) < |F - ip < C{\{F - /)wip + |(F - /)w 2 p + |(F - /)w 4 p) 
(11) = C{ 6 f + |(F - /)w 2 p + |(F - /)w 4 p) . 


By a simple geometric argument one finds 


( 12 ) |(F - /)w 2 p = ! + (! + 62 )^ - 2(1 + 62 ) cos ( 0 i 2 - f) , 

where 612 is the angle (measured anticlockwise) between wi and Fw 2 , which is 
determined by 


(13) 


cos 012 = 


(1 + ^ i )^ + (1 + , 52)^-(1 + ( 53 )^ 
2(1 + 5i)(1 + (52) 


and 


sin 012 > 0 , 


cf. [11, Proof of Lemma 2.2]. Notice that the condition sin 0 i 2 > 0 follows from the 
assumptions F G GL+(3) and Fw 2 • 62 > 0. 

Denote by 6 ij the acute angle formed by Fwi and Fwj and by 7744 that between 
Fw 4 and W4. Since 


(14) |(F - /)w 4 p = ! + (! + 154 )^ - 2(1 + (^ 4 ) COS 7744 , 


in order to express the right hand side of ( 11 ) in terms of the J^’s, we need to 
specialize cos 7744 in terms of the Si’s. Set 


Fw 4 


:= ( 01 , 02 , 03 ) 


|Fw 4 | 

and remark that, by assumption, 03 > 0. Thus 


(15) 


Fw 4 1 /3 /fi 

cos 7744 = 1 ^;^ ■ W 4 = 2“1 + 6 “2 + 3 

On the other hand oi and 02 are computed by solving 



(16) 

(17) 
where 

(18) 


Oi 


Fw 4 

- -- • Wi 

|Fw 4 | 


= cos 014 , 


oi COS 012 + 02 sin 012 = 


Fw 4 

|Fw 4 | 


Fw2 

|Fw 2 | 


COS 024 , 


COS 014 


(1 + <5i)^ + (1 + 64)'^ — (1 + Se)^ 

2(1 + ( 5 i)(l + ^4) 
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P5 



(19) 


cos 024 


(1 + 62)'^ + (1 + <54)^ — (1 + ^5)^ 
2(1 + (52)(1 + ^4) 


Taking into account (12)-(19), one can express the right hand side of (11) as a 
function f of d and see that /(O) = 0 and V/(0) = 0, which implies f{S) < C\5\^ 
for |(5| sufficiently small. For larger |i5| the inequality readily follows from (11)- 
(14). □ 


We will consider the octahedron O generated by the points Pi := (0, 0,0), P 2 := 
(1,0,0), P 3 := (0,1,0), P 4 := (1,1,0), P 5 := (^ ^ #), and Pg := (i,i,-^), 
see Figure 9. We call the triangulation determined by cutting O along the 
diagonal P 1 P 4 , further denotes the triangulation corresponding to P 2 P 3 , and 
the one corresponding to P 5 P 6 . 

Given a deformation u of the six vertices of O, denotes the piecewise affine 
extension of u corresponding to the triangulation i = 1,2,3. In the next 

lemma, Q denotes the interior of the (bounded) polyhedron determined by the 
images of the facets of O through any of the piecewise affine extensions u^'‘\ (Notice 
that the images of the facets do not depend on the chosen extension.) 

Lemma 3.3. One has that det > 0 a.e. in O and Q is convex if and only if 
det > 0 for every i = 1 , 2 ,3. 

Proof. Assume that det > 0 a.e. This implies that Q is connected, the diag¬ 
onal u{Pi)u{P 4 ) is contained in Q, and the outer normal vectors to the facets of Q 
point towards the outside of Q. Consider now the tetrahedra of the triangulation 
since the normals to the facets of Q point towards the outside of Q, it turns 
out that detVu^^) ^ 0 a.e. if and only if the diagonal 'u(P 2 )u(P 3 ) is contained in 
Q. The same holds for using the corresponding diagonal. On the other hand, 
an octahedron is convex if and only if all the three diagonals are contained in the 
inner part of the octahedron itself. □ 


The octahedron satisfies an estimate corresponding to the one of Lemma 3.2. 
Lemma 3.4. There exists C > 0 such that 
(20) dist^{Vu, SO{3)) <C “ 1)^ 

Pi,Pj NN 


a.e. in O , 
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for every u € such that u is piecewise affine with respect to the triangu¬ 

lation determined by cutting O along the diagonal PiPa, det Vm > 0 a.e. in O, and 
u{ 0 ) is convex. 

Proof. Let xc * = be the characteristic functions of the four tetrahedra 

Ti := P 1 P 2 P 4 P 5 , T 2 := PiP 2 ^ 4 ^ 6 i ^3 := P 1 P 3 P 4 P 5 , T 4 := P 1 P 3 P 4 PQ, respectively. 
Since Vu = J2i=iXiPi for some Fi G GL+(3), it suffices to prove (20) in each 
tetrahedron. Notice that Pi and P 4 are not nearest neighbours and therefore we 
cannot directly apply Lemma 3.2. On the other hand, the length of u{PiP 4 ), which 
is a common edge of the four deformed tetrahedra, can be expressed as a function of 
all the edges of u{0), the latter being a (possibly irregular) octahedron. Specifically, 
from the rigidity of convex octahedra, it follows that there exists a function / such 
that 

|u(PiP 4 )| = fih ,. •., I 12 ), 

where k, i = 1,..., 12, are the lengths of the twelve edges of u{0). In particular 
we set 

h := \uiPiP 2 )\ , h := \uiP 2 P 4 )\ , ^3 := 14^2^5)1 , h ■= |w(AP 5 )| , h ■= \uiP 4 P 5 )\ ■ 

The explicit formula of / is not important. Let Si := h — 1 for i = 1,..., 12, and 
So := |u(PiP 4 )| — \/2. We claim that / is differentiable at (1,..., 1). Then 

/(I + (5i,..., 1 + S 12 ) = \/2 + 0(|(5|), with (5 = ((5i,..., ^ 12 ), 

which yields, in combination with Lemma 3.2, the following inequality for Vu = Fi 
on the tetrahedron Ti: 

dist2(Fi,50(3)) < C (^(Z, - 1)2 + (|u(PiP4)| - 72 ) J = CY^SI 

5 2 

= CY,S‘Ip g(/(1 + ..., 1 + <5i2) - 72 ) < G|5|2 

for \ 6 \ < 1. On the other hand, by the triangle inequality, we have for \d\ > 1 

^0 < 2(f (1 + 5i,..., 1 + ^ 12 ) + 2) < Ailf + ll + 2 )< C'\S\^ . 

The inequality for the other Tfs is completely analogous. 

We are left to show that / is differentiable at (1,..., 1). To this end, we prove the 
existence and continuity of all its partial derivatives at (!,...,!). By a symmetry 
argument, it is enough to study the existence and continuity of d^f and 94 / (with 
reference to Figure 10a). 

We begin with dof. Let u{Pi) = Qi be such that h = 1 for each i ^ 3 and 
I 3 ^ 1. Since /^(l, 1, Z 3 ,1,..., 1) = 2 — 2 cosa, where a is the acute angle formed by 
QiQ 2 and Q 2 Q 4 , /(1,1, Z 3 ,1,..., 1) is a smooth function of a for 0 < a < tt. Next 
remark that the points Q 2 , Q 3 , Qs, Qe are coplanar; let p denote the plane containing 
them and LI be the orthogonal projection onto p (see Figure 10b). Considering the 
projections we see that 

|Q 2 n(Qi)| = |g3n(Qi)| = |Q5n(Qi)| = |g6n(Qi)| = cosf. 

Let 7 denote the acute angle formed by ^ 211 (^ 1 ) and Q 6 n(gi) (which is equal to 
that formed by ^ 611 (^ 1 ) and ^ 311 (^ 1 ) and that formed by ^ 311 (^ 1 ) and ^ 511 (^ 1 )). 
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Q2 


Qe 

(A) 


Qe 


Figure 10. (a) The image of O through a piece-wise affine map 
u such that h = 1 for each i ^ 3. (b) The projection of u{0) on 
the plane p, where O = n((5i) = n((54). 

Then, by the cosine formula 



( 21 ) 


and 




Remark that since a < tt, Z 3 is a smooth function of a. Moreover, since the derivative 

dh . 

—— is not zero at the point a = 7r/2, by the implicit function theorem it follows 
da _ 

that Z 3 is an invertible function of a in a neighbourhood of a = 7r/2, and its inverse 
is smooth in a neighbourhood of ^3 = 1. The differentiability of / with respect to 
I 3 at ( 1 ,..., 1 ) then follows from its smooth dependence on a. 

Finally, proving the existence and continuity of 94 / at (1,..., 1) is equivalent 
to proving that the length of Q 3 Q 2 is a smooth function of I 3 in a neighbourhood 
of (1,..., 1). The latter follows equivalently to the previous argument taking into 
account that IQ 3 Q 2 I = 2cosy- D 

Remark 3. Estimates (10) and (20) are crucial in the proof of the compactness of 
sequences of deformations with equibounded energy, as well as in the study of the 
T-limit and its scaling properties (see Theorem 4.1 and Proposition 1). Indeed, as 
already remarked, each of the lattices introduced in Section 1 defines a tessellation 
of the space into tetrahedra and octahedra. This allows us to deduce the following 
lower bounds on the energy 


(23) 

(24) 


dist^(VuE, S'0(3)) < C£^{ue) a.e. in 
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for each admissible deformation u^. Observe, in particular, that in the case of the 
diamond cubic lattice the above inequalities are obtained by first neglecting in the 
energy the bonds between atoms of the sublattice and then applying (10) and 
(20) on the tessellation of the space thus defined. 

Inequalities (23)-(24) imply, via the rigidity estimate (9), that Vue is locally 
close to SO{3) in and to ^S'0(3) in 


4. Dimension reduction and scaling properties of the F-limit 

In the present section we show the results in the three dimensional setting that 
were obtained in two dimensions in our previous paper [II]. The proofs of these 
results follow the lines of those in [11] by application of Lemmas 3.2 and 3.4. We 
will therefore omit further details of the proofs here. 

Given R G SO{3), p G (0,1], and fc G N, we define the minimum cost of a 
transition from an equilibrium in ^lk,oo with < 0 to an equilibrium in ^k,oo with 
> 0 as 

'y^{p,k,R) := inf {£^{v,p, k): M >0, v G Ap,oo{^k,oo) , 

Vu = I for xi G (—oo, —M) , Vv = ^R for xi G (M, +oo)} , 

where is defined in ( 8 ). 

Remark 4. In fact it can be proved that j^(p, k, R) does not depend on the rotation 
R (see [11, Proposition 2.4]). We therefore set 

7^(p,fc) ■=-f^{p,k,I) ■ 

The function 7 '^(p, k), which depends on the number of planes of atoms of the two 
lattices £]" and £+ contained in the domain £lk,oo, is in fact the relevant quantity 
that describes the system when e tends to zero. More precisely, our goal is to show 
that for k sufficiently large, there holds 

inf < 7^(1, fc), 

pG(O.l) 

i.e., the system displays dislocations. In order to prove this, we perform a dimension 
reduction with respect to the directions V 2 , V 3 . To this end, for each G Ap^e(S^ke) 
we define the rescaled deformation 

Ue{x) := Ue{AsX) , 

where is the matrix defined by 

{ AsYi = vi 

AeV2 = eV2 
AeV 3 = eV 3 . 

This yields a scaling of the domain £lke to VLk, which is independent of e. For fixed 
p G (0,1] and k G N we address the question of the F-convergence of the sequence 
of functionals {1^} defined by 
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where Ap^ei^k) is the corresponding set of admissible deformations, i.e. 

:= {■«£ G piecewise afhne, 

Vue constant on flk H (Aj^T) VT G Tp^J , 
det Vme > 0 a.e. in flk , 

Us{Aj^'P) is convex WV G Tp,e} ■ 


Theorem 4.1. Let {fte} C Ap^ei^k) be a sequence such that 

limSUpIe(Me) < C . 

E->0 + 

Then there exists a subsequence (not relabeled) such that 

VheA-^ ^ (SviM I d2 I da) weakly* in , 

where the functions u G R^), ^2,^3 G L°“(nfc;R^) are independent of V2 

and V 3 , i.e., dv^u = dvid2 = dv^d^ = 0, for i = 2,3. Moreover, 

~i , I ,x fco(5'0(3)) a.e.inVtf, 

{d,Md 2 \d^)Gl ' . A 

lco(^oO(3)) a.e. 

The sequence of functionals {X^} T-converges, as e ^ 0 “^, to the functional 

X{u) = ifuGA, 

I +c» otherwise, 

with respect to the weak* convergence in R^), where 

A := {u G R^): dv^u = dv^u = 0 a.e. in , 

| 9 viit| < 1 a.e. in Llf , I^vimI < - o,.e. in O^} . 


Remark 5. As in the two-dimensional case, the limit of a sequence of discrete 
deformations with equibounded energy does not depend on the triangulation chosen 
for the octahedra, see [11, Remark 3.4]. Moreover, the F-limit does not depend on 
the choice of the triangulation T'p^J in (7), since its formula only depends on the 
discrete values of the deformation, and not on its extension to the three-dimensional 
continuum. Similarly, it does not depend on the choice of the tessellation of . 


The next two results characterise the behaviour of the F-limit in the dislocation- 
free case and in the case when dislocations are present. 

Proposition 1 (Estimate in the defect-free case, p = 1). There exist Ci,C 2 > 0 
such that for every /c G N 

Cik^ <7^(l,fc) < C2k^. 

Proposition 2 (Estimate for p = X). There exist positive constants C[,C 2 such 
that for every k 

C[k'^ < 7^(A,fc) < 

Remark 6. The proof of Proposition 1 is a generalisation of its two-dimensional 
counterpart (see [11, Proposition 2.5]) given the rigidity estimates of the previous 
section. In contrast, the proof of Proposition 2 is straightforward: it follows by 
testing on the identical deformation v{x) = x and taking into account that each 
interfacial atom has a number of bonds that is uniformly bounded in k. 
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Remark 7. Propositions 1-2 in combination with Theorem 4.1 prove that dislo¬ 
cations are energetically preferred if the thickness of the nanowire modeled by k is 
sufficiently large. 
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